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ABSTRACT 

Context. Physical characterization of Trans-Neptunian objects, a primitive population of the outer solar system, may provide con- 
straints on their formation and evolution. 

Aims. The goal of this work is to characterize a set of 15 scattered disk (SDOs) and detached objects, in terms of their size, albedo, 
and thermal properties. 

Methods. Thermal flux measurements obtained with the Herschel-FACS instrument at 70, 100 and 160 yum, and whenever applicable, 
with Spitzer-MIFS at 24 and 70 yum, are modeled with radiometric techniques, in order to derive the objects' individual size, albedo 
and when possible beaming factor Error bars are obtained from a Monte-Carlo approach. We look for correlations between these and 
other physical and orbital parameters. 

Results. Diameters obtained for our sample range from 100 to 2400 km, and the geometric albedos (in V band) vary from 3.8 % to 
84.5 %. The unweighted mean V geometric albedo for the whole sample is 11.2 % (excluding Eris); 6.9 % for the SDOs, and 17.0 % 
for the detached objects (excluding Eris). We obtain new bulk densities for three binary systems: Ceto/Phorcys, Typhon/Echidna and 
Eris/Dysnomia. Apart from correlations clearly due to observational bias, we find significant correlations between albedo and diam- 
eter (more reflective objects being bigger), and between albedo, diameter and perihelion distance (brighter and bigger objects having 
larger perihelia). We discuss possible explanations for these correlations. 

Key words. Kuiper belt: individual: SDOs-detached objects - Infrared: planetary systems - Methods: observational - Techniques: 
photometric 



1. Introduction 

Since the detection of the first Trans-Neptunian Object (TNO) 20 years ago -1992 QBi- (Jewitt and Luu 1992) besides Pluto, about 
1550 TNOs/Centaurs have been discovered (as of October 2011). Diflferent dynamical classes, according to their orbital param- 
eters, have been identified within this distant icy population: classical objects with moderate eccentricity orbits; resonant objects 
trapped in mean motion resona nces with Neptune (Plutinos, in the 3:2 resonance with Neptune, represent the most populated reso- 
nant group, Jewitt & Luu 2000 ): scattered disk objects (SDOs) with high eccentricities and sometimes high inclinations, resulting 
from close encounters with Neptune; detached objects, whose perihelia are well outside Neptune's orbit CGladman et al. 2002 J , 

Send offprint requests to: P. Santos-Sanz: pablo.santos@obspm.fr 

* Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important 
participation from NASA. PACS: The Photodetector Array Camera and Spectrometer is one of HerscheVs instruments. 
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and which likely constitute a large population in the outer Kuiper Belt; centaurs, which reside between the orbits of Jupiter and 
Neptune, represent TNOs scattered towards the inner solar system, and so are usually considered as Kuiper Belt Objects (KBOs). 
Although it is convenient to classify objects by "populations", the limits between some of these populations (i.e. between SDOs, 
centaurs, high inclination/high eccentricity classicals, detached, etc..) are often "blurred" and may ultimately rely on definition 
(.Gladman et al. 2008 ). The present paper deals with 15 SDO or detached objects, for which we first specify the adopted definition. 

The terminology "scattered disk" was originally used for TNOs with orbits characterized by both a high eccentricity and a 
perihelion close to Neptune. It is believed that the majority of SDOs have unstabl e orbits due to clo se encounters with Neptune, but 
some of them could have survived since the origin of the solar system (Dunca n & Levison 1997b . Alternatively, the origin of the 
scattered disk could be related to other eff'ects, including sweeping resonances that move the objects to this region jMalhotra 1995K 
passing star (or stars) fM orbidelli & Levison 20041 and references therein). Several approach es to the classification of objects a s 
SDOs numerically integrate their orbits backward in time over dozens of millions of years (Elli ot et al. 2 005 , Gladmai Tet al 2008]) . 
Here we adopt the Gladman et al. (2008) classification, which defines SDOs as those objects, independently of their origin, that are 
currently scattering actively off' Neptune. SDOs are characterized by the rapid variation of semi-major axis in the orbit integrations. 
Specifically, Gladman at al. (2008) adopt a criterion similar to Morbidelli et al. (2004): the object is classified as an SDO if its 
semimajor-axis presents a variation > 1.5 AU during a 10-m.y. backward integration. The exact value used (1-2 AU variation) 
makes little diff'erence, as SDOs undergo large semimajor-axis changes in short times. Then, SDOs in this nomenclature might 
be called "scattering" objects instead of "scattered" objects. Following this classification SDOs exist over a large semi-major axis 
range, not only restricted to a > 50 AU; rather the SDO population extends down to the Centaur region (beginning at a = 30 AU). 
At the other edge the inner Oort cloud would begin at very large semimajor-axes (a > 2000 AU, Gladm an et al. 2 008 ). 

Detached objects dDelsanti & Jew itt 2006 ) are also known as "extended scattered disk objects" (ESDOs), though it is preferable 
to use the term "detached" because not all these objects have necessarily been placed in this region through scattering. Detached 
objects have orbits with perihelion distances unaff'ected ("decoupled") from the giant planets, in particular Neptune. This population 
may include inner Oort cloud bodies or transitional objects between SDOs and inner Oort cloud objects. The limits in parameter 
space are unkno wn, and are expan ding with the new discoveries of TNOs at larger and larger semimajor axes -e.g. a= 965 AU 
for 2006 SQ372- (Gladman et al. 2002 ). Elliot et al. (2005) adopted the Tisserand parameter (T) relative to Neptune to define the 
detached population. A problem with this definition is the large number of high inclination TNOs within the classical population. 
In the Gladman et al. (2008) scheme, numerical integrations separate SDOs from detached TNOs. A remaining issue is to decide 
where the detached population ends towards low eccentricities; Gladman et al. (2008) adopt a lower bound of ^ = 0.24. Possible pro- 
cesses that placed the detached objects on these orbits are currently being investigated (e.g. lGomes et al. 20081 Kenyon et al. 2008 
IDuncan et al72Q()8l etc). 

Determining KBO sizes and albedos is not a straightforward task. Measurements from direct imaging are available for only a 
handful of th em (Bro wn & Trujillo 2004[IBrown et al. 20061) . Most recently, sizes have become available from stellar occultations 
dEUiot et al. 2010t |Sicardy et al. 201l| )rObtaining size/albedo values for a larger sample is best achieved from thermal radiometry, 
whereby thermal and optical observations are combined within a thermal model. After the pioneering measurements on two objects 
from ISO dThomas et al. 2 000 ) and ground-based observations on a few others (see |Stansberry et al. 2008|and references therein), 
Spitzer provided the first large such dataset, with about 60 Centaurs/KBOs measured at 24 and/or 70 //m dCruikshank et al. 20051 
Grundy et al. 2005 & 2007, Stansberry et al. 2006 & 2008, IBrucker et al. 20091) . This approach (multi- wavelength photometry) 
also constrains the th ermal regime of the object, phenomenologically described by the "beaming factor" (Jones & M orrison 19741 
[Lebofsky & SpenceT 1989 , Harris 199 8 ). The launch of Herschel Space Observatory in 2009 off'ered a new tool for this purpose, 
providing photo met ry over 6 bands fro m 70-500 yum. Early results on a handful of objects were presented in IMiiller et al. 20101 
ILim et al. 20 10' and' Lellouch et al. 2010 Beyond their intrinsic scientific value, albedos are needed to properly interpret reffectance 
spectra. Furthermore, correlations between size, albedo, color, composition and orbital parameters may provide diagnostics on the 
dynamical, collisional and physical history of KBOs. The size distribution of TNOs is also a diagnostic relating to that history, but 
in the absence of measured sizes the distribution has typically been estimated from the luminosity function and an assumed albedo, 
an uncertain step given that the average albedos seem to be diff'erent for the diff'erent dynamical populations ( [Grundy et al. 2005 



IBrucker et al. 20091) . Instead, and to constrain Kuiper belt formation/evolution models it is desirable to directly determine the size 
distribution of the diff'erent populations, and infer their individual total masses. Large bodies retain information about the angular 
momentum of the presolar nebula and their size distribution is dominated by accretion. On the other hand, small bodies -in the 50 
to 100 km range- are the result of collisional evolution ( P etit et al. 2 008i 

A large program (370 h of telescope time), named "TNOs are Cool: A Survey of the Transneptunian Region" 
dMtiller et al. 20091) , has been awarded on Herschel to measure the thermal ffux of about 140 KBOs/Centaurs. We present here 
results on 15 SDO and detached objects. These observations lead to new sizes/albedos for 9 objects and improved results on 6 
others. Results on 19 classical and 18 Plutinos are presented in companion papers by Vilenius et al. (2012) and Mommert et al. 
(2012), respectively. 



2. Observations 

The Herschel Space Observatory d Pilbratt et al. 2010l) is a 3.5 m telescope launched on 14 May 2009, covering the far-infrared range 
by means of three focal plane instruments. The "TNOs are Cool" program consists of photometry of -140 objects, mostly with the 
Photodetector Array Camera and Spectrometer (PACS, Poglitsch et al. 2010 ). 

The PACS instrument has 3 photometric bands centered at 70, 100 and 160 jdm. These bands (termed 'blue', 'green' and 'red', 
respectively) cover the 60-85 yum, 85-130 jum, 130-210 /im ranges. With two bolometer arrays, PACS can observe a field in two 
bands at a time: blue/red or green/red. The red detector is a 32 x 16 pixel array, while the blue/green is 64 x 32 pixels. The 3 
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bands image the same field of view of 3.50' x 1.75'. Projected onto the sky, pixel sizes are 6.4'' x 6.4" (red band) and 3.2" x 3.2" 
(blue/green bands), respectively. 

The observations discussed here were performed using the scan map mode. This mode was selected due to its better overall 
performance compared to the point-source (chop-nod) mode (Miill er et al. 20101 . Scan maps are performed by slewing the space- 
craft at a constant speed along parallel lines, or "legs". Our observations were typically performed at medium speed (20"/^). We 
use 10 scan legs in each map, separated by 4". The length of each leg is 3.0'. Each such scan-map is repeated two to five times, 
depending on the expected brightness of the objects (see below). Within a single scan map, either the blue-red or the green-red 
channel combination is selected. For any band selection, two individual maps are made using two scanning directions with respect 
to the detector array: 70° and 110°, allowing to fit the scan pattern on the rectangular array in an optimum way. 

The "TNOs are cool" target list was constructed by predicting the thermal fluxes in the PACS bands for the known TNOs using 
the "hybrid standard thermal model" (STM - see Sect.|4]). The flux predictions were based on absolute visual magnitudes (H) from 
the Minor Planet Center (MPC), an assumed visual geometric albedo of 8% (based on previous studies using Spitzer), and geometric 
circumstances appropriate to the expected observing epoch (c. 2010). For objects with known albedos (for this study those were 
Eris: 70%, Typhon: 5%, 1996 TL66: 4%, and 2002 PN34: 5%) we used those values rather than the 8% assumption. Objects with the 
largest predicted thermal fluxes were selected as potential targets; that list was further refined by rejecting objects with predicted 
positional uncertainties < 5". The objects presented in this paper are a subset of the targets selected in the above manner. Once the 
targets were selected, integration times were chosen to provide signal-to-noise ratio (SNR) values of ~ 10 in the blue and green 
channels. Integration times for individual observations were tailored by varying the repetition factor for the scan maps described 
earlier. Note that the values of H from the MPC were used for planning purposes only, and can diff'er significantly (by up to 0.5 
mag) from the values ultimately used in our analysis. 

It is important to note that the target selection process introduces a bias towards brighter (low H values) objects. A consequence 
of the assumed and generally fixed albedo is that the object sample favors large objects at large heliocentric distances (r/^). In 
addition, the object sample is aff'ected by a discovery bias that: (i) favors high-albedo objects at high r/^, and (ii) introduces a 
correlation between the perihelion (q) and the heliocentric distance, because objects are most easily discovered near perihelion 
and haven't moved much since discovery. These biases and correlations must be kept in mind when studying correlations between 
size, albedo and other physical and orbital parameters (see Sect. [5]). As for the present SDOs/Detached sample, objects are listed in 
Table [T] including relevant orbital parameters, absolute magnitudes in V-band, rotational periods and amplitudes, taxonomic color 
class (when known), existence or not of satellites, and dynamical classification following Gladman et al. (2008). Observational 
circumstances are described in more details in Table |2l 

Each target was visited twice, with identical observations repeated in both visits. During the first visit, we grouped 4 maps to 
be observed in sequence, corresponding to the two scan directions and the two channel combinations. The observation geometry is 
such that on the final combined map, the central ~ 50" x 50" region has high coverage. Photometric accuracy is usually limited 
by instrumental noise in the blue channel and by background confusion noise in the red, while either of the two sources of noise 
may dominate in the greeiQ The timing of the obser vations, i.e. the selection of the visibility window, was optimized to ensure the 
lowest far-infrared confusion noise circumstances (K iss et al. 2005 ), in such a way that the estimated SNR due to confusion noise 
has its maximum in the green channel. The timing of the second visit was calculated such that the target moved 30-50" between 
the two visits. This ensures that the target position during the second visit is still within the high-coverage area of the map from the 
first visit. At the same time, this enables us to determine a nd s ubtract the background in the individual maps, so as to obtain a clean 
combined image where we extract the final flux (see Sect. |3.1| ). 

3. Data reduction and photometry 

In this section, we first describe the reduction process applied to the PACS data in order to obtain what we call "single" and 
"combined" maps, and then the photometric techniques used to determine the target fluxes and associated uncertainties in each 
band. Further details of the data reduction are given in Appendix A. 

3.1. Data reduction 

The data reduction of the PACS maps was performed within the Herschel Interactive Processing Environment (HIP^ version 
6.0.2044) by means of adapted standard HIRE scripts. Since our targets move slowly (a few "/hour), and the total observation time 
per scanning direction in one visit is only 10-25 min, we do not correct for the apparent motion. 

After the application of the HIRE scripts we obtain one "single" map per visit, filter, and scan direction (i.e. in total 8 maps per 
object in the red, and 4 maps in the blue/green). The final re-sampled pixel scale of the single maps are 1.1 "/pixel, 1.4"/pixel, and 
2.1"/pixel for the blue (70 yum), green (100 jim) and red (160 jim) channels respectively. 

We use these single maps to generate final "combined" maps on which the photometry will be performed. To generate the 
combined maps, we combine the two visits to the target to determine the background map and subtract it from each single map. 
Finally we co-add all the background- subtracted maps in the co-moving frame of the target. Similar methods were used to subtract 
background sources in Spitzer TNO observations, and are described in Stansberry et al. (2008) and Brucker et al. (2009). An 
alternate, simpler method, to obtain final maps is to co-add directly the original single maps (i.e. not background-corrected) in the 

^ PACS AOT Release Note: PACS Photometer Point/Compact Source Mode, 2010, RICC-ME-TN-036, Version 2.0, custodian T. Muller, |http : | 
[TZherschel . esac . esa . int/twiki/bin/view/Public/PacsCalibrationWeb 

^ HIRE is a joint development by the Herschel Science Ground Segment Consortium, consisting of ESA, the NASA Herschel Science Center, 
and the HIFI, PACS and SPIRE consortia members, see: http://herscheLesac.esa.int/DpHipeContributors. shtml. 
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Table 1. Orbital parameters, absolute magnitudes, B-R colors, photometric variation, taxonomy, dynamical classification and mul- 
tiplicity of the observed objects 



Object 


a 


q 


i 


e 


Hy 


B-R 


P 


Avclr 


Taxon. 


Class. 




[AU] 


[AU] 


[deg] 




[mag] 


[mag] 


[h] 


[mag] 






1996 TL66 


84.5 


35.0 


23.9 


0.59 


5.39+0. 12^'^'^^ 


1.11+0.07^''?''-'^ 


12.0^ 


<0.12, <0.06^'^ 


BB 


SDO 


2001 FPi85 


212.0 


34.3 


30.8 


0.84 


6.39+0.07^'^ 


1.40+0.06^'^ 




<0.06^ 


IR 


SDO 


2002 PN34 


31.2 


13.4 


16.6 


0.57 


8.66+0.03^'^ 


1.28+0.02^ 


8.45/10.22^^ 


0.18+0.04^ 


BB-BR 


SDO 


2002 XU93 


66.5 


21.0 


77.9 


0.69 


8.11+0.10"'^ 


1.20+0.02^ 








SDO 


2007 ORio 


67.1 


33.6 


30.7 


0.50 


1.96+0.16^ 










SDO 


2007 RWio 


30.4 


21.2 


36.0 


0.30 


6.39+0.6F 










SDO 


(2003 FX128) Ceto*'^ 


99.7 


17.8 


22.3 


0.82 


6.54+0.06^' 


1.42+0.04^ 


4.43« 


0.13+0.02" 




SDO 


(2002 CR46) Typhon*^ 


37.6 


17.5 


2.4 


0.54 


7.72+0.04^'^'^'^'^'^'"' 


1.29+0.07"^'^ 


9.67/19.34^'« 


0.06+O.OP'^ 


BR 


SDO 


1999 KR16 


48.7 


33.9 


24.9 


0.30 


5.37+0.08"^ 


1.87+0.07^^'"'" 


5.93/11.86" 


0.18+0.04" 


RR 


DO 


2003 FY128 


49.2 


37.0 


11.8 


0.25 


5.09+0.09^ 


1.65+0.02^ 


8.54^ 


0.15+O.OF 


BR 


DO 


2005 QU182 


114.0 


37.0 


14.0 


0.68 


3.80+0.32^ 










DO 


2005 TB190 


76.5 


46.2 


26.4 


0.40 


4.40+0.11^ 


1.54±0.03^ 








DO 


2007 OCio 


50.0 


35.5 


21.7 


0.29 


5.43+0.10^'' 










DO 


2007 UKi26*^ 


74.0 


37.7 


23.4 


0.49 


3.69+0.10^ 










DO 


(2003 UB313) Eris*^ 


68.1 


38.6 


43.8 


0.43 


-1.12+0.03^^'^ 


1.21+0.09^'^ 


13.69/28.08^ 


<0.1,0.01^'^ 


BB 


DO 



32.13/25.92^'^ 



* Indicates that the object is a known binary/multiple system. Hy [mag]: average visual magnitude obtained from papers referenced below. 
Numbers in bold indicate Hr magnitudes instead of Hy. See text for details about the absolute magnitude estimation. B-R [mag] colors (when 
available). P [h] single or double-peaked rotational period. Ahir [mag]: lightcurve amplitude. Taxon.: taxonomic color class (Perna et al. 2010 
and references ther ein). Class.: dynamical cl assification fo llowing Gladman et al. 2008. SDO = sc attered disc object, DO - detached ob ject. 
References: Ro manishin & Tegler 2005] IDoressoundir am et al.2005 , Jewitt & Luu 2001] IRabinowitz et al. 20071 "^[Sheppard 2 010, 
Derived from MFC V-data; Derived from MF C R-data; Benecchi et al. 2009 . Perna & D otto 2011, priv. c o mm.; ^^[Grundy et al. 2008, 
IPerna etal. 20 10 '"^ Peixinho et al. 2004, Shep pard & Jew itt 2002, Belskay a^etal. 2003[ ^^[Sheppard 2007[ [Tegler & Romanishin 1998) 
IJewi tt& Luu 1998 Barucci et al. 1999, Tegler et al. 2003, ^^Trujillo & Brown 2002[ ^ IBoehnhardt et al. 2002[ Tegler 2011, priv. comm.; 
^^lOrtiz et al. 2006, Thirouin et al. 2010 , Sheppar d & Jewitt 2003[ ^^ lOrtiz et al. 20031 ^^ IDuff^ard et al. 20081 ^^ IRoeetal. 20081 
[Grundy et al. 2007 1 ^ ^Brown & Schaller 2007i : ^^ iNoll etal. 20091 



Table 2. Individual observational circumstances 



Object OBSIDs Dur. [min] Mid-time rn [AU] A[AU] a[dQg] 



(15874) 1996 TL66 


1342190953 


-0956 


39 


7 


22-Feb-2010 10:57:34 


35 


7245 


35 


8103 


1 


59 


(15874) 1996 TUe 


1342191029 


-1032 


39 


7 


23-Feb-2010 12:09:47 


35 


7250 


35 


8289 


1 


59 


(82158) 2001 FPi85 


1342211422 


-1425 


97 


4 


23-Dec-2010 08:41:10 


34 


9353 


35 


2646 


1 


53 


(82158) 2001 FPi85 


1342211528 


-1531 


97 


4 


23-Dec-2009 21:36:49 


34 


9356 


35 


2562 


1 


53 


(73480) 2002 PN34 


1342213067 


-3070 


41 





18-Jan-2011 14:18:48 


17 


5092 


17 


8769 


2 


99 


(73480) 2002 PN34 


1342213089 


-3092 


41 





18-Jan-2011 21:31:02 


17 


5098 


17 


8821 


2 


98 


(127546) 2002 XU93 


1342204211 


-4214 


59 


8 


09-Sep-2010 14:56:03 


21 


1063 


21 


4432 


2 


58 


(127546) 2002 XU93 


1342204240 


-4243 


59 


8 


09-Sep-2010 21:19:45 


21 


1064 


21 


4399 


2 


58 


2007 ORio 


1342220081 


-0084 


78 


6 


07-May-2011 02:16:22 


86 


3285 


86 


5992 





65 


2007 ORio 


1342220272 


-0275 


78 


6 


09-May-2011 23:59:21 


86 


3305 


86 


5534 





66 


2007 RWio 


1342213219 


-3222 


59 


8 


24-Jan-2010 05:43:55 


27 


4552 


27 


8385 


1 


90 


2007 RWio 


1342213270 


-3273 


59 


8 


24-Jan-2011 21:05:07 


27 


4558 


27 


8495 


1 


89 


(65489) Ceto 


1342202877 


-2880 


59 


8 


ll-Aug-2010 20:15:24 


31 


6508 


31 


8552 


1 


81 


(65489) Ceto 


1342202910 


-2913 


59 


8 


12-Aug-2010 11:00:59 


31 


6524 


31 


8666 


1 


80 


(42355) Typhon 


1342210596 


-0599 


59 


8 


3O-N0V-2OIO 21:13:47 


18 


2160 


18 


5731 


2 


90 


(42355) Typhon 


1342210624 


-0627 


59 


8 


Ol-Dec-2010 03:50:07 


18 


2162 


18 


5690 


2 


91 


(40314) 1999 KR16 


1342212814 


-2817 


97 


4 


17-Jan-2011 20:35:59 


35 


7593 


36 


0633 


1 


51 


(40314) 1999 KR16 


1342213071 


-3074 


97 


4 


18-Jan-2011 15:51:40 


35 


7589 


36 


0497 


1 


51 


(120132) 2003 FY128 


1342212770 


-2773 


59 


8 


16-Jan-2011 19:01:01 


38 


8406 


38 


6763 


1 


45 


(120132) 2003 FY128 


1342213107 


-3110 


59 


8 


19-Jan-2011 05:03:50 


38 


8416 


38 


6356 


1 


44 


2005 QU182 


1342212619 


-2622 


41 





14-Jan-2011 14:46:12 


48 


8994 


49 


1349 


1 


13 


2005 QU182 


1342212696 


-2699 


41 





15-Jan-2011 14:01:49 


48 


9008 


49 


1522 


1 


12 


(145480) 2005 TB190 


1342221729 


-1732 


78 


6 


27-May-2011 03:17:15 


46 


3134 


46 


5135 


1 


24 


(145480) 2005 TB190 


1342221782 


-1785 


78 


6 


28-May-2011 13:15:54 


46 


3132 


46 


4900 


1 


25 


2007 OCio 


1342206671 


-6674 


59 


8 


16-Oct-2010 21:52:04 


35 


5187 


35 


1842 


1 


54 


2007 OCio 


1342206698 


-6701 


59 


8 


18-Oct-2010 11:48:52 


35 


5189 


35 


2090 


1 


55 


(229762) 2007 UK126 


1342202277 


-2280 


41 





08-Aug-2010 18:09:36 


44 


9448 


45 


1702 


1 


27 


(229762) 2007 UK126 


1342202324 


-2327 


41 





09-Aug-2010 13:05:14 


44 


9440 


45 


1572 


1 


27 


(136199) Eris 


1342199487 


-9490 


59 


8 


30-Jun-2010 21:01:08 


96 


6522 


96 


8807 





59 


(136199) Eris 


1342199753 


-9756 


59 


8 


03-Jul-2010 19:45:58 


96 


6518 


96 


8326 





60 



OBSIDs: Herschel internal observation IDs. Dur. [min]: total duration of the 4 OBSIDs observations for the red band. For the blue and green 
channels the time is half of this value. Mid-time: Mean date and UT time of the observation. Vh [AU]: heliocentric distance (AU) at mid-time. A 

[AU]: distance objQci-Herschel (AU) at mid-time. Qf[deg]: phase angle -degrees. 
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target frame. This method is obviously less optimal than the previous one in terms of SNR, but provides a useful test to demonstrate 
that the background subtraction is not introducing any spurious effects. A more complete and detailed description about the data 
reduction process can be found in Appendix A, and will be published in Kiss et al. (in prep.-di). 



3.2. Photometry 

We use IRAF/Daophoj^ to perform standard synthetic-aperture photometry on the final sky-subtracted, co-added maps for our 
targets. The target is typically obvious and located very near the center of the maps. For very faint targets and/or ones with some 
residual structure from the background subtraction (as can happen in the green and red channels), the ephemeris position of the 
target is used to position the photometric aperture. 

Once the target is identified, we measure the flux at the photocenter position for aperture radii ranging from 1 to 15 pixels. We 
perform an aperture correction technique dHowell 1989b for each aperture radius using the encircled energy fraction for a point- 
source for the PACS instrumenj^ We construct an aperture-corrected curve-of-growth using these fluxes (red line of Figure A.4). 
Given the optimized reduction technique above, the fluxes used to construct these aperture-corrected curves-of-growth need not be 
corrected for background contribution. 

After inspection of the curves-of-growth, the optimum synthetic aperture to measure the target flux is selected, and is typically 
~ 1.0-1.25 times the point spread function (PSF) FWHM in radius (PSF FWHM in radius is "^l" 111" I XT' in blue/green/red bands 
respectively). Usually, the optimum aperture radius lies in the "plateau" zone of the curves-of-growth (see Figure [O) . In cases in 
which there are other sources or artifacts very close to the target, we use a smaller aperture radius to avoid contamination of the 
target flux. 

Uncertainties on the flux measurements are estimated by means of a Monte-Carlo technique, in which 200 artificial sources are 
implanted on the actual background- subtracted final maps. These artificial sources have the structure of the PSF for the specific 
(blue, green, or red) channel. Sources are randomly implanted on a square region of 50'' x 50'' around the target photocenter, with 
an exclusion zone corresponding to a circular region around the target photocenter with a radius - 2xPSF FWHM. We measure and 
aperture-correct the fluxes of these artificial sources using an aperture radius of 5 re-sampled pixels, which is the median optimum 
aperture radius. Uncertainties are computed as the standard deviation of these 200 fluxes, and finally multiplied by a factor V2. The 
latter factor is due to the fact that given the way the final maps are generated, the remaining background is measured only once in 
the immediate vicinity of the real target but two times in the rest of the image (i.e. for the artificial sources). Note that except for the 
background subtraction, all the above reduction and photometric steps were validated on standard stars with known flux. 

The fluxes and associated uncertainties are finally color-corrected for each band. Flux densities within the PACS photometric 
system are defined as those that a source with a flat spectrum would have at the PACS reference wavelengths of 70, 100 and 160 
\im (Poglitsch et al. 2010). Converting the in-band flux to monochromatic flux densities therefore depends on the spectral energy 
distribution of the source and the transmission profile for each PACS band. Color-corrections for each filter are tabulated as a 
function of color temperatur^ but are usually within a few percent except at the very low end {J < 30 K) of the temperature 
range. For the color temperature, we used an approximation of the disk-averaged temperature, i.e. Tcoior - T^^/2^/'^, where T^^, the 
subsolar temperature (see equation ([2]) ), was obtained from a first run of the thermal model described below, using at this step the 
in-band (i.e. color-uncorrected) fluxes. T color values for our sample range from 25 K for Eris to 77 K for 2002 PN34. The mean color 
correction factors for aU our sample are: 0.991+0.028 (blue band), 0.989+0.006 (green band), and 1.011+0.017 (red band). The 
impact of these factors is smaller than the uncertainties of our measurements. 

The absolute flux calibration of PACS is based on standard stars and large main belt asteroids and has uncertainties of 3% for 
the blue/green channel, and 5% for the red channe|^ The final flux uncertainties we use in our modeling are the root- sum- square 
combination of these absolute calibration uncertainties and the uncertainties we obtain from the photometric measurements for 
individual targets. A summary of the final, color-corrected and absolute-calibrated, fluxes of our SDOs/Detached sample is given in 
Tables 



3.3. Spitzer-MIPS observations 

About 102 TNOs/Centaurs have been observed with the 0.85 m Spitzer Space Telescope dWerner et al. 20041) using the Multiband 
Imaging Photometer (MIPS) instrument (Rieke et al. 2004t . The MIPS 24 jdm channel consist of a 128 x 128 pixel detector and 
the 70 yum channel of a 32 x 32 pixel detector. For these 2 channels the telescope-limited resolution is 6" and 18" respectively. 
Fluxes for 60 objects have been published in Stansberry et al. (2008) and Brucker et al. (2009) but others are unpublished. Eight 
targets of our sample have available (published or unpublished) Spitzer-MIFS flux density measurements at 24 and/or at 70 yum (see 
Table |4]). Unpublished data have been processed in a similar way as in Stansberry et al. (2008) and Bru cker et al. (2009), using the 



calibration ( Gor don et al. 20071 |Eng elbracht et al. 2007 ) and color-corrections ( Stan sberry et al. 2007| ) described in the references. 



As in Stansberry et al. (2008), we adopt absolute calibration uncertainties of 3% and 6% for the 24 and 70 fim observations (50% 



^ Image Reduction and Analysis Facility (IRAF), is a collection of software written at the National Optical Astronomy Observatory (NOAO) 
geared towards the reduction of astronomical images in pixel array form. 



Miiller et al. 2011: PACS Photometer -Point-Source Flux Calibration, PICC-ME-TN-037, Version 1.0; Retrieved November 23, 2011; |http: 
|77herschel . esac . esa . int/twiki/pub/Public/PacsCalibrationWeb/pacs_bolo_f luxcal_report_vl . pd£ 

^ Miiller et al. 2011, PACS Photometer Passbands and Colour Correction Factors for Various Source SEDs, PICC-ME-TN-038, Version 1.0; 
Retrieved November 23, 201 1; http : //herschel . esac . esa . int/twiki/pub/Public/PacsCalibrationWeb/cc_report_vl . pdf 

6 Miiller et al. 2011: PACS Photometer -Point-Source Flux Calibration, PICC-ME-TN-037, Version 1.0; Retrieved November 23, 2011; |httpr| 
[//herschel . esac . esa . int/twiki/pub/Public/PacsCalibrationWeb/pacs_bolo_f luxcal_report_vl . pdf | 
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Table 3. Color corrected fluxes at 70, 100 and 160 [im for the SDOs/Detached sample 



Object 


Flux^^ Tm Jvl 


FluXinnrmJvl 


Fluxi AnrmJvl 


Class. 


(15874) 1996 TL66 


9.6+1.2 


9.7+1.1 


9.2+2.2 


SDO 


(82158) 2001 FPi85 


9.2+0.8 


11.1 + 1.1 


7.0+1.9 


SDO 


(73480) 2002 PN34 


11.4+1.2 


11.9+1.4 


7.6+6.1 


SDO 


(127546) 2002 XU93 


16.5+0.9 


15.0+1.6 


6.2+2.1 


SDO 


2007 ORio 


3.7+1.1 


4.8+1.4 


6.7+2.3 


SDO 


2007 RWio 


13.2+0.9 


11.3+1.3 


10.3+2.1 


SDO 


(65489) Ceto 


12.3+0.7 


10.5+1.0 


6.2+4.8 


SDO 


(42355) Typhon 


28.5+1.2 


21.4+1.2 


9.1+2.4 


SDO 


(40314) 1999 KR16 


5.7+0.7 


3.5+1.0 


4.6+2.2 


DO 


(120132) 2003 FY128 


14.4+0.9 


14.0+1.2 


13.0+2.2 


DO 


2005 QU182 


4.5+0.9 


2.5+1.1 


8.4+3.0 


DO 


(145480) 2005 TB190 


5.8+0.9 


8.6+1.4 


5.4+1.8 


DO 


2007 OCio 


8.3+0.9 


7.8+1.2 


7.4+1.8 


DO 


(229762) 2007 UK126 


13.7+1.1 


12.3+1.3 


7.2+2.2 


DO 


(136199) Eris 


2.0+0.6 


3.7+1.0 


5.6+1.6 


DO 



Fluxvo, Fluxioo? Fluxi6oi Color-corrected fluxes (mJy) at 70, 100 and 160 jdm. Class., dynamical classification following Gladman et al. 2008. 

SDO = scattered disc object, DO = detached object. 

Table 4. Spitzer-M\PS color corrected fluxes at 23.68 and 71.42 jdm available for our SDOs/Detached sample. 



Object 


Mid-time 


rJAU] 


A [AU] 


Qf[deg] 


FIUX24 [mJy] 


Flux7o[mJy] 


(15874) 1996 TUe 


2008-Feb-17 04:27:11 


35.4370 


34.9378 


1.4 


0.324+0.020' 


10.12+1.40" 


(82158) 2001 FRi85 


2004-Jun-20 22:51:52 


34.2594 


34.0927 


1.7 


0.395+0.096^ 


Contaminated flux 


(73480) 2002 PN34 


2006-Jul-16 11:00:56 


14.6074 


14.1524 


3.6 


10.144+0.1172 


30.34+1.592 


(127546) 2002 XU93 


2OO6-N0V-O4 20:09:23 


21.0896 


20.8251 


2.7 


2.431+0.0523 


14.12+1.193 


(65489) Ceto 


2006-Jul-19 03: 19:41 


27.9924 


27.6753 


2.0 


1.419+0.0472 


14.6+2.6^ 


(42355) Typhon 


2008-Jun-26 22:18:45 


17.6784 


17.5110 


3.3 


4.974+0.1253 


27.46+3.432 


(120132) 2003 FY128 


2007-Jul-14 07:30:00 


38.3556 


38.1281 


1.5 


0.427+0.05 P 


19.38+3.163 


(136199) Eris 


2005-Aug-24 22:13:07 


96.9049 


96.4083 


0.5 


< 0.0051 


2.7+O.71 



Mid-time: Mean date and UT time of the observation, r [AU]: heliocentric distance (AU) for the observation date. A [AU]: distance 
oh]Qci-Spitzer (AU) for this date. Qf[deg]: phase an gle -degrees. FIUX24, F luxyp: Spitzer-MW^ color-corrected flux den sities (mJy) at 23.68 fim 
and 71.42 jim. Upper limits are Icr. References: |Stansberry et al. 2008[ 2^ Revised data from |Stansberry et al. 2008 1 3^ Previously unpublished 
data. 



larger than the uncertainties derived for observations of stellar calibrators). The larger uncertainties allow for efl'ects from the 
background- subtraction technique, the faintness o f the targets, and uncertainties in the color corrections. The background subtraction 
techniques used for our Herschel data (see Sect. |3.1| ) are in fact derived from techniques originally developed for the MIPS data 
reductions. The previously unpublished fluxes presented here are based on new reductions of the MIPS data, utilizing updated 
ephemeris positions for the targets. This allows for more precise masking of the target when generating the image of the sky 
background, and for more precise placement of the photometric aperture. In most cases the new fluxes are very similar to the 
previously published values for any given target, but in a few cases significant improvements in the measured flux and SNR were 
achieved. 

As discussed below, the combination of Spitzer-M\PS and Herschel-FACS fluxes significantly improves results from the thermal 
modeling. The Spitzer 70 pm band is similar to the Herschel 70 pm channel, providing a consistency check on data from the two 
observatories. The inclusion of the Spitzer 24 pm flux, when available, strongly improves the ability to constrain the surface temper- 
ature distribution, i.e. in general, allows us to determine the beaming factor (see Sect.|4]). Published, re-analyzed, and unpublished 
Spitzer-MWS flux density measurements available for our sample are given in Table |4] Note that the eflTective wavelengths for these 
fluxes are 23.68 and 71.42 pm. 

4. Thermal model 

4.1. General formulation 

The combination of thermal-infrared and optical data makes it possible to infer an objects' diameter (Z)), geometric albedo {py) and 
constrain its thermal properties. This requires the use of a thermal model to calculate the disk-integrated thermal emission. Here we 
use the so-called hybrid standard thermal model (hybrid STM, Miiller et al. 2010, Stansberry et al. 2008). Except for the fact that 
the hybrid STM assumes zero phase angle (which is an excellent approximation for TNOs), it is identical to the Near-Earth Asteroid 
Thermal Model (NEATM, Harris 1998 ) . NEATM was originally developed for Near Earth Objects (NEOs) but is applicable to all 
atmosphereless bodies and is being used in companion papers by Mommert et al. (2012) and Vilenius et al. (2012). In this model, 
the temperature distribution follows instantaneous equilibrium of a smooth spherical surface with solar input, but is modified by a 
factor ri (beaming factor), which represents empirically the combined effects of thermal inertia and surface roughness. As shown by 
Harris (1998) and Stansberry et al. (2008), multi- wavelength measurements can be used to constrain 7/ if the data are of suflftcient 
quality and encompass both the short-wavelength and the peak and/or long-wavelength portions of the spectral energy distribution 
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(SED) of the target. In this case the data give a direct measurement of the color-temperature of the emission, which is sensitive to the 
value of T]. Our preferred approach was to simultaneously model Herschel and Spitzer measurements (of course accounting for their 
different observer-centric, A, and heliocentric, r/^, distances) and fit an r] factor to the combined dataset. We call this the ''flo ating-/;" 
approach. When Spitzer data were unavailable (or when this method led to implausible y; values; see IMommert et al. 2012 ) we 



instead used a "fixed-z/" approach, with 7/ = 1.20 ± 0.35, the mean value derived by |Stansberry et al. 2008 from Spitzer observations 
of TNOs. 

The model therefore fits the Herschel (3 data points) or Herschel -h Spitzer (5 data points) fluxes in terms of two {D and py) 
or three (Z), py, and if) paramet ers. The model also requires the choice on a phase integral that we assumed to be related to py 
through q = 0.336 • py -h 0.479 dBrucker et al. 20091) . Diameter and geometric albedo are further related through: 

D = 2' 10^©/^ • 10"^/^ ^^P^ ' lAU/km (1) 

where H is the absolute magnitude (nominally in V band), and V© = -26.76 + 0.02 jBessell et al. 19981) the V magnitude of the Sun. 
A grey emissivity (e= 0.9) was assumed when calculating the local temperatures and monochromatic fluxes. With these definitions 
the subsolar temperature is given by: 

'Fq • (1 -py-q)^^^"^ 



(2) 



where ¥q is the solar flux at 1 AU. The solution for the D, py (and rj) parameters are determined from a minimization of the 
deviation of the model from the observed fluxes, accounting for their individual 1-cr error bars. Note that when only upper limits 
were available, they were treated as non-detections (i.e. zero flux) with a l-cr uncertainty as the upper limit. In the case of multiple 
systems D is the area-equivalent diameter, i.e., that of the sphere with the same total surface area, because none of the binaries are 
resolved in our data. 



4.2. H magnitudes 

Equation ([T]) above outlines the role of absolute magnitude (H) in the thermal modeling. Determining accurate values of H, prefer- 
ably in V band, is therefore an essential task for this project. We describe hereafter our adopted approaches. 

In the best situation, photometrically calibrated observations are available at a variety of phase angles. In this case, correcting 
the observed V magnitudes for the heliocentric (r/^) and geocentric (A) distances provides my(l, 1) = V - 5 - log(r • A). A linear fit 
of the my (1, 1) magnitudes vs phase angle a (m(l, 1) = Hy -\- ft - a), taking into consideration measurement errors, then provides 
the phase coefficient JS (mag/deg) and the absolute magnitude Hy, and their associated error bars. Whenever data are unsufficient 
in V band but adequate in R band, we similarly and alternatively infer Hr, which by means of the thermal model will yield the red 
geometric albedo (pr). 

In the case when not enough data in either V or R are available, or if the above approach fails due to poor quality data, we 
adopt a mean /3 = 0.10 ± 0.04 mag/deg which is the average value (excluding Pluto) from jBelskaya et al. 2008, We then derive H 
magnitudes for each individual photometric measurement and use the weighted average result. 

When photometrically-calibrated magnitudes are not available in the literature, we use astrometric observations compiled by 
MFC, as those include rough estimates of the optical magnitudes (most often in R band). The above approaches are then applied 
to this data, bearing in mind that i ndividual V or R magnitudes reported in MFC data may have large (and systematically biased) 
uncertainties up to 0.5 magnitudes ([Romanishin & Tegler 2005K 



When full lightcurve information is available (i.e. the amplitude is known and the lightcurve can be temporally phased to the 
Herschel measurements), the H magnitude is corrected accordingly. When the lightcurve amplitude (or an upper limit) is known, 
but the period is unknown or phasing is impossible, 88% of half peak-to-peak amplitude (i.e. Icr = 68% of the values of a sinusoid 
are within this range) is quadratically added to the uncertainty on H. Finally, if there is no lightcurve information, we assume peak- 
to-peak amplitude less than 0.2 mag, as based on Duffard et al. (2009), -70% of the objects have a peak-to-peak amplitude less 
than this. In this case we quadratically add 0.09 mag (i.e. 88% of 0.2/2) to the H uncertainty. Absolute magnitudes quoted in Table 
[TJhave been obtained as just described. 

Note finally that in case we have to use R magnitudes (i.e. for 2007 ORio, 2007 RWio and 1999 KRig) modeling provides the R 
geometric albedo (pr) instead of V albedo (py). Pr can be converted to py through: 



Pv = Pr 



. iQ[((V-R)e-(V-RUj))/2.5] 



where (V - R)q= 0.36 mag is the Sun's V-R, and (V - R)obj that of the object. We assume (V - R)obj= 0.55+0.11, the mean value 
from Hainaut & Delsanti (2002) for 66 SDOs/Detached objects, except for 2007 ORio (see discussion for this object in Sect.|5.3|). 



4.3. Uncertainties 

Within the framework of the thermal model, uncertainties on the fitted (D, py, rj) parameters were obtained using a Monte Carlo 
approach (see Mueller et al. 201 1) in which 1000 synthetic datasets were randomly generated using the uncertainties in the measured 
fluxes. Note that the uncertainty in the H magnitude (and the +0.35 uncertainty in rj, in the case of the "fixed-//" approach) was 
also accounted for in the Monte-Carlo approach. Naturally, we eliminated those Monte-Carlo simulations which generated negative 
fluxes (or negative //). The diameter, albedo (and rj) of each of the 1000 synthetic objects were determined, and their distributions 
were used to determine the error bars on these parameters. As outlined in Mueller et al. (201 1), these distributions are generally not 
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Fig. 1. Radiometric fixed-// fits for the objects observed only with Herschel-V ACS . 



Gaussian (especially the albedo distribution). Therefore, we adopted the median of the Monte-Carlo results as the nominal value, 
and asymmetric error bars to include the 68.2 % of the results. 
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Fig. 2. Radiometric hybrid-STM fits for the objects observed with Spitzer-MIFS and Herschel-V ACS. The best fit is plotted for the 
Herschel (r/^, A) conditions (solid line) and for the Spitzer conditions (dashed line). 



Figure [T] shows best fits for all objects in our sample using only the Herschel-FACS measurements. All these fits are "fixed-77". 
For the 8 objects that also have Spitzer fluxes, Figure [2] shows best fits ("floating-//", except for Eris, which is not detected at 24 jim) 
for the combined Herschel and Spitzer measurements. In that Figure, the same best fit model is plotted twice, once (solid line) for 
the Herschel (rh, A) conditions and once for the Spitzer (dashed line) conditions. In general, this shows that the model appropriately 
accounts for rather diff'erent Herschel vs. Spitzer fluxes associated to difl'erent rh and A (e.g. Ceto, 2001 FPigs). Many of the fits 
shown in Figure [T] and [2] are satisfactory, in that the model fits all data points within or close to the l-cr error bars. However, in 
a number of cases, the fits are obviously bad (e.g. 1999 KR16, 2005 QU182, Typhon); in those cases, the data points are mutually 
inconsistent and the model is no more than a mere compromise between them. Given that the PACS green vs blue observations (and 
the Spitzer vs Herschel data) are not taken simultaneously, such behaviour could in principle result from rotational variability of the 
objects. An alternate physical explanation would be that the fluxes are afl'ected by spectral features, e.g. depressing the 100-yum flux 
with respect to the 70-yum and 160-yum values, or more generally that our simple STM-like approach does not capture the complexity 
of the thermal emission. However, the discrepancies are often large and these explanations may not be quantitatively plausible. In 
those cases, we conclude that the error bars we have determined on the fluxes are probably underestimates of the true error bars, 
though we are unable to track down the reasons for that. 

In the situation where the formal measurement errors are too small (or similarly that the model is not completely adequate), 
the above described Monte-Carlo method will underestimate the uncertainties on the solution parameters. To handle this problem, 
we adopted a "rescaled error bar" approach described in detail in Appendix B (Sect. |B.1| ). Essentially, it consists of "inflating" 
the measurement error bars by a suitable factor before running the Monte-Carlo simulation. This approach was applied to the six 
objects for which the radiometric fit is poor (see Figs. [T] and [2]): 1999 KR16, 2002 PN34, 2002 XU93, 2005 QU182, 2007 UK126, and 
Typhon, with rescaling factors varying from 1.45 to 1.99. Fluxes for the remaining nine objects are satisfactorily fit within the 1-cr 
flux uncertainties and did not require any adjustment of the measurement errors. We finally stress that the diameter, albedo, and 
beaming factor values and their uncertainties reported in the next section are obtained in the framework of the hybrid STM, which 
may have its own limitations (see Harris 2006 for a discussion in the case of NEOs). More physical models (e.g. TPM) may lead to 
somewhat diff'erent values and error bars, as illustrated by Miiller et al. (2010) in the case of a few objects observed in the Science 
Demonstration Phase (SDP) of the "TNOs are cool" programme. 
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5. Results and discussion 

5.1. General results 

Results of the model fits in terms of the objects diameter, albedo and beaming factor are gathered in Table 5. When Spitzer-MIFS 
measurements are available (8 objects), the model was run twice, once considering the PACS fluxes only, and once combining 
PACS and MIPS. As illustrated in Table 5, these two fitting options give consistent solutions within error bars. The PACS+MIPS 
solution is always adopted as the preferred one (outlined in bold face in Table [5|, as it (i) provides reduced error bars (ii) permits 
the determination of the beaming factor. 

Table [5] includes first size/albedo determination for 9 objects (4 SDO and 5 detached objects). Six of them are based on PACS- 
only measurements (2007 ORio, 2007 RWio, 1999 KRi6, 2005 QU182, 2007 OCio, and 2007 UK126), and the remaining three make 
use of PACS plus unpubHshed-MIPS data (2001 FPigs, 2002 XU93, and 2003 FY128). PubHshed sizes and albedos exist for the other 
6 objects of our sample (1996 TLgg, 2002 PN34, Ceto, Typhon, 2005 TB190 and Eris) but we here present fluxes for these objects 
at wavelengths not observed before (i.e. 100 jim and 160 jum), and the combination with earlier measurements leads to improved 
estimations of sizes, albedos, and beaming factors (see Table[5]). 

The V geometric albedos for our targets vary from 3.8 % (2002 XU93) to 84.5 % (Eris), with an unweighted mean value of 
11.2+7.6 % excluding Eris. If weighted by the relative errors (i.e. l/(cr/pv)^), the mean albedo is 6.4+6.3 % (for cr we used the 
mean of upward and downward uncertainties on py). This error bar on the mean albedo includes the dispersion of the albedo 
measurements and their individual RMS errors. It is dominated by the dispersion between individual measurements, so it has a 
limited significance. This variation of geometric albedos can be seen in Figure [3]-left. For SDOs, there is a lack of albedos>20 % 
and a clustering in the [0-5 %] range. In contrast, detached objects in our sample all have albedos> 5%, with a peak in the [10-20 
%] range. The unweighted/weighted mean is 6.9/5.2 % for the SDOs and 17.0/12.3 % for detached objects excluding Eris. 

Diameters vary from 112 km (2002 PN34) to 2454 km (Eris), with a remarkably uniform distribution (Figure |3] right). However, 
the distribution of diameters is seemingly diff'erent between SDOs and detached objects. SDO diameters are smaller than 400 km, 
except for an outlier at 1280 km (2007 ORio), while all detached objects are larger than 200 km. This apparent tendency of detached 
objects to be larger than SDOs is, however, probably an eff'ect of the discovery bias (see discussion in Sect. [2]), as detached objects 
are typically discovered at larger heliocentric distances than SDOs. 

Except for Eris, where the "floating 77" fit fails, the combination of Spitzer and Herschel data allows us to determine the beaming 
factor for 7 objects (6 SDO and 1 detached -see Table [5]-). The (weighted by absolute errors) mean-// obtained from these results is 
1.14+0.15. This is fully consistent with the mean 1.20+0.35 //-value derived from the Spitzer dataset (Stansberry et al. 2008) and 
justifies the use of this mean value for those cases where 7/ cannot be determined. Beaming factors only barely larger than ~1 (as 
opposed to 7/ » 1) indicate that KBOs in general and our sample objects in particular are close to the slow rotator regime (i.e. 
STM), which is a remarkable result given the very low temperatures prevailing on these objects and their usually short rotation 
periods. Except in the unlikely situation where they are all seen pole-on, this implies very low thermal inertia and/or very large 
surface roughness. Lellouch et al. (2011) determine the thermal inertia of Charon's surface to be F = 10-20 J m"^s"^/^K"^ (MKS) 
and find this to be equivalent to 7/ = 1.21 - 1.40. Yet Charon is a slow rotator (6.4 day) and a rough scaling indicates that the same 
ri range, applied e.g. to Haumea or Orcus, would be representative of even much smaller thermal inertias (F = 0.5-3 MKS). Hence 
a general conclusion is that unless they are very rough, KBO surfaces must have very low thermal conductivity, suggesting a high 
degree of surface porosity. Similar thermal inertias and t/'s for KBOs were obtained in Miiller et al. (2010). 
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Fig. 3. Left: geometric albedo histograms for SDOs/Detached, SDOs and detached objects. Right: diameter histograms for 
SDOs/Detached, SDOs and detached objects. 



5.2. Comparison with earlier results 

Six objects of our sample (1996 TLge, 2002 PN34, Ceto, Typhon, 2005 TB190 and Eris) have previously published sizes/albedos. 
We here briefly compare results for the first five, and discuss Eris in more detail in next subsection. Note that 2007 UK126 has been 
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reported as binary system dNoU et al. 2009b , but there is no mass estimation for this object, for this reason so we do not discuss it 
here. 

(15874) 1996 TL66. Our preferred solution for this SDO (D = 339+20 km, pv = n.O^J. \ %, ?] = 1.15+^-^^) is at odds with the 



-1.5 'I ^'^^-0.05^ 

fll6 Vm ri.. - ^ <+2.0 , 



Stansberry et al. 2008| results (D = 575+}}^ km, pv = 3.5+^-^ %, 7]=1.H± 0.3) based on Spitze r-MIFS fluxes at 24 a nd 70 jim. Our 



reprocessing of these data leads to very diff'erent fluxes at 24 jim (0.32 mJy vs 0.38 mJy in jStansberry et al. 2008| ) and especially 
at 70 yum (10.1 mJy vs 22 mJy). Fitting these revised Spitzer fluxes alone gives a best fit {D= 328 km, py= 11.8 %, 7]= 1.12) in 
excellent agreement with our preferred solution. 

(73480) 2002 PN34. Although the 70 jdm fluxes from Spitzer and Herschel are very difl'erent (Figure [2|), this can be ascribed to 
the difl'erent observing distances. Our preferred results are: D = 112 + 7 km, py = 4.9 + 0.6 %, 77 = 1.02^ ^^, in a good agreement 
with [Stansberry et al. 2(jQ8| results: (D = 120 + 10 km, pv = 4.3^^-^ %, rj = l.lO^^ j^). This object is the smallest one in our 
SDOs/Detached sample. 

(65489) Ceto/Phorcys is a binary SDO, previously observed with Spitzer-MIFS ( [Grundy et al. 2007| [Stansberry et al. 2008| ). 
Model fits indicate that the optimum model (Z) = 281 + 11 km, pv = 5.6 + 0.6 %, j] = 1.04 + 0.05) matches the 24 jdm flux but 
somewhat overestimates the Spitzer 70 yum fl ux, indicative of a cold er object than the Spitzer data would suggest. It is therefore 
not surprising that the Spitzer-only solution of [Stansberry et al. 2008 has a lower 77 (0.86^^^^), lower D (230^ |g km) and higher pv 
(7.7^|;| %). 

Grundy et al. (2007) estimate a mass of (5.41+0.42)xl0^^ kg for the Ceto/Phorcys system. Using an equivalent diameter of 
218^22 ^^^y obtain a bulk density of 1.37^^32 g cm"-^. We derive a smaller bulk density of 0.64^q|3 g cm"^ due to our larger 
equivalent diameter (Z) = 281 km which gives Dceto = 223+10 km, and Dphorcys = 171 + 10 km, assuming the same albedo for the 
two components and Am = 0.58+0.03 mag). This low bulk density is compatible with a porous surface, which is also supported by 
the very low t] = 1.04 value. 

(42355) Typhon/Echidna is a binary SDO. Our preferred results are Z) = 185+7km,;7y = 4.4+0.3%, 77 = 1.48+0.07, consistent 
with the Stansberry et al. (2008) solution obtained by applying a fixed 77 = 1.2+0.35 model to the 24 jim flux (D = 174^|g km, 
Pv = 5.1^Qg%). In contrast, these results are inconsistent with the preliminary Herschel results (D = 138 + 9 km and pv = 8.0 +1.0 
%), based on SDP data (Miiller et al. 2010 ). The latter observations were performed in chop-nod mode (less sensitive than scan-map 
for point sources) and led to lower 70 and 100 jim fluxes and no detection at 160 jim. 

Grundy et al. (2008) obtain a mass of (9.49+0. 52)x 10^^ kg for this binary system. Using an equivalent diameter of 179^|g km, 
they obtain a bulk density of 0.44^^^^ g cm"^ . We derive a nominally even smaller (but compatible within error bars) bulk density 
of 0.36^QQ^ g cm"^ from a slightly larger equivalent diameter (D = 185 km, i.e. Dxyphon = 162 +7 km, and DEchidna =89 + 6 km, 
assuming same albedo and Am = 1.30+0.06 mag). Again, this extremely low bulk density points to a porous body. 

(145480) 2005 TB190. This detached object was not observed by Spitzer-MWS, but was already measured in chop-nod mode by 
Herschel. Our new results, using a fixed-77 mode (J) = 464 + 62 km, pv = 14.8^3-^%) are reasonably consistent (to within 1-cr) with 
the preliminary values pubHshed by Miiller et al. (2010) (D= 375 + 45 km, pv = 19.0 + 5.0%). 



5.3. Two large objects, 2007 ORio and Eris 

2007 ORio. Its Hr magnitude of ~2 indicates that this very distant object (currently at 86.3 AU) must be one of the largest TNOs, 
unless it exhibits a very high albedo. Brown et al. (201 1) recently reported that its spectrum is markedly red in the near-infrared and 
shows features due to water ice. This is an unsual combination as most KBOs with water ice absorptions have neutral colors, but 
this characteristic is shared by Quaoar. Although their data ar e unsufficient to show a detection of methane ice. Brown et al. (201 1) 
used the analogy with Quaoar, where methane ice is present ( Schaller & Brown 2007a ), to suggest that the red color of 2007 ORio 
is due to irradiation of small amounts of methane ice. They further assumed an albedo equal to that of Quaoar (nominally 0.18 for 
an adopted Quaoar radius of 960 km, but allowing for 50 % uncertainties), and showed, based on the Schaller and Brown (2007b) 
sublimation model, that 2007 ORio is in a similar volatile retention regime as Quaoar, with the main volatile species (N2, CO, CH4) 
being marginally stable over the solar system age. 

2007 ORio is clearly detected in the three FAGS channels. In the lack of V measurements, we assumed Hr = 1.96 + 0.16 
magnitudes based on the fit of eight good quality R-measurements from MFC data (see Table [T]). A fixed-// fit is satisfactory (Figure 
[TJ), providing D = 1280 + 210 km, pr = 18.5^^ 2%. This puts the object in the Charon, Makemake, Haumea, etc... size class, i.e. 
a plausible c andidate for dwarf pla net status. To convert pr to pv we use equation ^ from Sect. 42_ but use Quaoar's V-R value 
= 0.67+0.02 JFornasier et al. 20041) . This yields pv = 13.9^3 ^%, ^i^^iy consistent with the value assumed by Brown et al. (2011). 
2007 ORio thus indeed appears as a sibling to Quaoar, probably of even larger size. 

(136199) Eris. Although Eris was observed by Spitzer, it was not detected at 24 yum (with a 0.005 mJy upper limit), preventing 
us from determining its beaming factor. Assuming therefore 77= 1.20 + 0.35, we obtain a best fit solution for D = 2454 + 117 km, 
Pv = 84.5 + 8.8 %. Using the mass determined by Brown et al. (2007) -recomputed in Grundy et al. (201 1)- for the Eris/Dysnomia 
system ((1.688+0.035)xlO^^ kg ) we obtain a bulk density of 2.23^^33 g cm"^, assuming the same albedo for Eris and Dysnomia 
(DEris= 2434 + 117 km, and Doysnomia^ 316 + 23 km). It is probably more realistic to assume that Dysnomia's surface is much 
darker than Eris'. Assuming for definiteness a 5 times lower albedo, we obtain Z)£n^=2356 + 1 17 km and DDysnomia= 685 + 50 km, 
yielding a global density o f 2.40^|j;3^ g cm~-^, similar to the Brown et al. (2007) value (2.3 + 0.3 g cm"^). 

Near-IR spectroscopy ([Tegler et al. 2010', Merlin et al. 20091 [Brown et al. 2005) indicates that Eris' surface composition is very 



similar to Fluto's, with the presence of methane ice along with another dominant ice, presumably nitrogen, suggestive of a high 
albedo. Frior to our m easurements, Eris' si ze and albedo have been determined by several groups and methods. The first measure- 
ment was achieved bv lBertoldi et al. 20061 who used IRAM thermal observations at 1200 jum and obtained D = 3000^^^^ km and 
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Table 5. Radiometric fit results for the objects sample. These are the results obtained by the application of an hybrid- STM to the 
fluxes obtained by Herschel-FACS (with Spitzer-MIFS data when are available). Preferred results for each target are in bold 



Object 



Data 



D [km] 



Pv [%] 



(15874) 1996 TL66 
(15874) 1996 TLge 
(82158) 2001 FPi85 
(82158) 2001 FPi85 
(73480) 2002 PN34 
(73480) 2002 PN34 
(127546) 2002 XU93 
(127546) 2002 XU93 
2007 ORio 
2007 RWio 
(65489) Ceto 
(65489) Ceto 
(42355) Typhon 
(42355) Typhon 
(40314) 1999 KR16 
(120132) 2003 FY128 
(120132) 2003 FY128 
2005 QU182 
(145480) 2005 TB190 
2007 OCio 
(229762) 2007 UK126 
(136199) Eris 
(136199) Eris 



PACS 
PACSAIIPS 

PACS 
PACS/MIPS 

PACS 
PACSAIIPS 

PACS 
PACSAIIPS 

PACS 

PACS 

PACS 
PACSAIIPS 

PACS 
PACSAIIPS 

PACS 

PACS 
PACSAIIPS 

PACS 

PACS 

PACS 

PACS 

PACS 
PACS/MIPS 



341+40 
339+20 
329+42 

332:^1 
110+12 

112+7 
173+20 

164+9 
1280+210 
247+30 
300+39 
281+11 
175+17 

185+7 

254+37 

479+55 

460+21 
416+73 
464+62 
309+37 
599+77 
2420+;?o 
2454+117 



11.0^ 
11.0^^ 

4.6+0.7 

5 2+1-3 
4.9+0.6 

o c+0.9 
•^•^-0.6 

3.8+0.4 

18.5+J-^* 
8.3!6-8* 

A Q+1-6 

5.6+0.6 

4 9+11 

4.4+0.3 

'■^-1.5 

7.9+1.0 
32.8 
14.8 

12.7+^-^ 
16.7:^1 

84.5+8.8 



Class. 



16.0 
10.9 
+5.1 



1.20+0.35* 

1 -|c+0.08 
-^•-^''-0.05 

1.20+0.35* 

1.20+0.35* 
1 02+0-07 

A-"^_0.09 

1.20+0.35* 

1 -t^+O.OS 
A.l^_0.08 

1.20+0.35* 
1.20+0.35* 
1.20+0.35* 
1.04+0.05 
1.20+0.35* 
1.48+0.07 
1.20+0.35* 
1.20+0.35* 
1.07+0.08 
1.20+0.35* 
1.20+0.35* 
1.20+0.35* 
1.20+0.35* 
1.20+0.35* 
1.20+0.35* 



SDO 
SDO 
SDO 
SDO 
SDO 
SDO 
SDO 
SDO 
SDO 
SDO 
SDO 
SDO 
SDO 
SDO 
DO 
DO 
DO 
DO 
DO 
DO 
DO 
DO 
DO 



Object: SDO or detached object observed with Herschel-FACS (and Spitzer-MIFS when available). Data: PACS, model results using only fluxes 

fi-om Herschel-FACS. PACS/MIPS, model results using fluxes from Herschel-FACS and Spitzer-MIFS. D [km]: Diameter and uncertainty, in 
kilometres, obtained from the thermal model. Pv [%]: geometric albedo and uncertainty, in V-band, obtained from the thermal model. * indicates 
geometric albedos derived from Hr magnitudes instead of Hy, (i.e, [%]). 77: Parameter from the thermal model application. * objects for which 
there are not Spitzer-MIFS fluxes available, in this case we assume a mean t] value to perform the radiometric fits. Class., dynamical 
classification following Gladman et al. 2008. SDO = scattered disc object, DO = detached object. 



Pv = 60+ 10 %. Subsequentlv JBrown et al. 2006lestimated D = 240 0+100 km, = 86 + 7 % from direct imaging with the Hubble 
Space Telescope. Based on the Spitzer 70 jim flux, |Stansberry et al7 2QQ8 obtained D = 2657^209 ~ ^^-^-in Spitzer- 

MIFS and hybrid-STM fit. Within uncertainties, our results agree with all these values except with those of .Bert oldi et al. 20061 
which overestimate the size and underestimate the albedo. Most recently, a stellar occultation by E ris on November 6th 2010 pro- 
vided a very accurate determination of Eris' diameter and albedo: D = 2326 ±12 km, pv = 96^^ % ( [Sicardy et al. 201 1) , assuming 
a spherical shape. (Note however that allowing for elliptical shape, Sicardy et al. find a broader range of solutions, 68.3 % (l-cr) of 
which have an efiTective diameter of 2330 ±180 km). Eris' diameter is thus indistinguishable from Pluto's given error bars on both, 
supporting the idea that Eris is a Pluto twin, with an even brighter surface, presumably due to a collapsed atmosphere at aphelion. 

At face value and given error bars, our result on Eris' diameter agrees nicely with that of Sicardy et al. (2011), especially if 
we assume that Dysnomia is intrinsically much darker than Eris. Similarly, our density for the system is in line with theirs (2.52 
± 0.05 g cm"^). However, the consistency is actually not complete. Sicardy et al. (2011) reanalyzed the Spitzer 70 pm and IRAM 
1200 pm fluxes in the light of their determined diameter and albedo. They considered both a standard thermal model (STM) and an 
isothermal latitude model (ILM), introducing an "y/"(= tjr) parameter describing exclusively surface roughness efifects (and ranging 
from t]r = l.O -no roughness- to ijr = 0.7 -large roughness-). In this framework, the only free parameter of their thermal model 
was the phase integral q. They found that the thermal measurements implied q = 0.49-0.66 for the STM case (fully consistent with 
the value of Saturn's brightest satellites), but an implausible q = 0-0.24 for the ILM. On this basis, the STM was strongly favored, 
implying either a pole-on orientation or a negligible thermal inertia. The inconsistency with our analysis resides in the fact that 
we assumed r]= 1.20 ± 0.35 and that our pv solution implies q = 0.75, according to the assumed Bruc ker et al. (2009) expr ession. 
Another way of seeing this is to note that if we now assume the two limiting STM cases considered by [Sicardy e t al. 20TT](i.e. q= 
0.49, T]R = 1.0, and q= 0.66, tjr = 0.7) and re-fit the PACS + MIPS data, we obtain best fit diameters of 2034 km and 2139 km 
respectively, and V- albedos larger than 100% in the two cases, inconsistent with the occultation results. Ultimately, the discrepancy 
can be traced to the fact that Sicardy et al. (2011) included the IRAM 1200 pm flux in their model, which as they noted tends to 
indicate warmer temperatures than the Spitzer 70 pm flux, while on the contrary our best fit model tends to underestimate the latter. 
The resolution of this problem will require a simultaneous fits of all datasets, considering in particular (i) non- spherical shapes 
(ii) separate thermal models for Eris and Dysnomia and (iii) more elaborate thermophysical models in which the thermal inertia, 
surface roughness and spin orientation are free parameters. This study is deferred to a subsequent paper (Kiss et al. in prep.-h). Note 
finally that the interpretation of thermal observations of a volatile-covered body such as Eris is expected to be even more complex 
than presented here. An additional complication is that volatile transport will result in surface variegation, i.e. multiple surface 
albedos. This is demonstrated in the case of Pluto, which does show large albedo spatial variations and whose thermal data require 
three-terrain model (ILellouch et al. 20111 IMonmiert et al. 20121 . Another example is Makemake, the visual appearance of which is 
unknown but for which a two-terrain model is required to fit the radiometric data (ILim et al. 20101 . 
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5.4. Search for correlations 
5.4.1 . Description of the method 



Using Spearman-rank correlation coefficient ( [Spearman 19Q4| ), also known as Spearman-p, we have examined our results for corre- 



lations between albedo (py), beaming factor (77), diameter (D), absolute magnitude (//), color (B-R), perihelion distance (q), orbital 
inclination (/), semi-major axis (a), and heliocentric distance at the time of observation (r/^). 

Spearman-p is distribution-free and less sensitive to outliers than most other methods but, like others, treats data-points as 
'exact' and does not take into account their possible error bars. Any variations in the measured data points within their error bars, 
may change the correlation coefficient, however. Furthermore, each correlation coefficient has its own confidence interval, which 
depends on the number of data points and on the magnitude of the correlation. We use procedures described in Appendix B (Sect. 
|B.2[ ) and in Peixinho et al (2004) to deal with these eff'ects to obtain Spearman-p values that take into account the error bars. 

We call correlations "strong" when |p| > 0.6. Efron & Tibshirani (1993) analyze the evidence for correlation with the follow 
criteria (where CL is the confidence level defined in Appendix B): CL ^ 99% (p < 0.010): very strong evidence of correlation; and 
CL > 99.7% (p < 0.003), the 'standard' 3cr significance: clear correlation. We restrict ourselves to the cases with CL > 99.5% (p < 
0.005), i.e. at least 2.8(T confidence. 




Fig. 4. Plots illustrating the three more interesting correlations with CL ^ 99.5% (^ 2.8cr) in our results. From left to right: albedo 
v^- diameter, albedo perihelion distance, and diameter v^- perihelion distance (Eris is out of the albedo scale in the first two panels, 
but its albedo is indicated). 



5.4.2. Correlations analysis 

Based on the previously described methods, the strongest and/or more interesting correlations are shown in Table|6] These are given 
for the whole sample (including or excluding Eris), and for SDOs and detached objects separately. Before discussing them, we stress 
that results must be taken with care because of small number statistics. In general, more significant correlations are found for SDOs 
than for detached objects. 

"Clear correlations" (> 99.7%) in Table |6] are seen between D and H (negative correlation, 3. Oct), py and H (negative, 4.2(t), 
D and rh (positive, 4.9(t), pv and rh (positive, 3.8(t). The first two anticorrelations are a simple consequence of equation (1), with 
no physical implication. The positive correlations between D and rh, and between pv and rh, are most probably biases due to object 
selection criterion and discovery bias, as discussed in Sect. [2] 

Additional "clear" or "very strong evidence" (> 99%) for positive correlations are found between D and q (3 Act) and p^ and q 
(2.9(t) (Figure]?]). There is no a priori relationship between diameter, albedo and perihelion distance, hence these correlations could 
plausibly be real, as discussed below. However, they could also reflect biases due to the correlation between q and rh (at 3.9cr for 
our sample, see Table [6]), i.e. the discovery bias discussed in Sect. [2] On the other hand, it is noteworthy that no high albedo objects 
(e.g. Py >10 %) are seen at small heliocentric or perihelion distance (<30 AU), which should not be a bias. Note that when SDOs 
and detached objects are treated separately, the D ys q correlation remains significant at 2cr only, while the py and q correlation 
essentially vanishes. Correlations also decrease (marginally) if Eris is excluded from the sample. 

Based on Spitzer measurements of about 50 objects, Stansberry et al. (2008) noted the apparent positive correlation between py 
and q (3.5cr in their case). While recognizing the above discovery bias, they proposed it to be real, off'ering two explanations. First, 
objects near the Sun tend to experience higher temperatures, hence may be subject to increased sublimation of their intrinsically 
bright ices. Conversely, increased UV-photolysis and solar wind radiolysis associated with the low heliocentric/perihelion distances 
may preferentially darken the surfaces of these objects. It seems more difficult to find a physical explanation for an hypothetical 
correlation between diameter and perihelion distance, especially because the current orbits of the SDOs/Detached objects tell us 
little if anything about the formation region. This apparent correlation is most probably a bias, as discussed above. Note also that 
unlike Stansberry et al. (2008) we do not find any (even tentative) correlation between py and the semi-major axis a (p= 0.21, 0.8cr, 
see Table |6]). 

We find a very strong correlation between albedo and diameter (pv vs. D) when SDOs are included, in particular for the whole 
sample at 2.9cr significance (see Figure]?]). This result was obtained by Stansberry et al. (2008) for KBOs (defined there as any 
object but Centaurs) but was not robust against small changes in the classification. This correlation implies a trend for larger objects 
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having higher albedos. An explanation could invoke the fact that larger objects can retain bright ices (released either by collision 
or intrinsic activity) more easily than small objects do. Eris follows this tendency, although its extremely high albedo puts it as an 
outlier in the plot of Figure]?] An important aspect here is that there are no small objects (<300 km) with high albedo (>10 %). This 
cannot be a bias, as such objects would have been targetted in our sample, yet poorly detected. Note that if objects with albedos 
>10 % are to be associated with the presence of surface ices, then such ices should be only moderately volatile (i.e. CH4 or less), 
as volatile retention models dSchaller & Brown 2QQ7b[|Brown et al. 201 II) show that only the largest (> 1000- 1500 km) bodies can 
retain the most volatile ices such as N2 or CO over the solar system age. 

Finally, we do not detect any correlation between diameter/albedo and color, or inclination, nor between the beaming factor 77 
with any parameter (though the sample is limited to 7 objects in the latter case). 

In summary, in spite of the relatively small sample covered in this paper, we find three possible correlations (i.e. pv vs. D, pv 
vs. q, and D vs. q), of which at least the first one -large objects tend to have brighter surfaces than small objects- is real and may 
be physically explained. More subtle correlations may be found, but with lower significance, and larger datasets will be required for 
confirmation. 

6. Summary 

We have determined albedos/sizes of 8 SDOs and 7 detached objects using Herschel Space Observatory-PACS fluxes measurements 
at 70, 100 and 160 fim and Spitzer-MWS fluxes at 23.68 and 71.42 iim when available. The main results, which include nine first 
determinations, are summarized hereafter: 

- Diameters range from 100 to 2400 km. The overall diameter distribution in our sample is remarkably uniform, but apparently 
diff'erent for SDOs vs detached objects. This diff'erence might however just reflect a sample/discovery bias, as our sample favors 
large objects at large heliocentric distances, and detached objects are typically discovered at larger distances than SDOs. 

- V-geometric albedos range from 3.8 to 84.5 %. Unweighted mean albedo for the whole sample is 11.2 % and 6.9/17.0 % for 
the SDOs/Detached (excluding Eris) respectively. There is a lack of albedos > 20% with a peak in the [0-5%] range for SDOs. 
Detached objects have albedos > 5 %, with a peak in the [10-20 %] range. 

- We re-determine bulk densities of three binary systems, based on our estimated sizes and published masses: Ceto/Phorcys 
(0.64+^1^ g cm"^), Typhon/Echidna (0.36+^-^^ g cm"^), and Eris/Dysnomia (2.40+^ ^6 g cm"^). 

- Beaming factor (7/), determined from thermal fits of 7 objects, range from 1.02 to 1.48, with a weighted mean obtained of 
1.14 + 0.15, consistent with the mean value derived from Spitzer, 1.20 ± 0.35 (Stansberry et al. 2008). These low 77 values are 
compatible with low thermal inertias, presumably due to high porosity of regolith-like surfaces. 

- We obtain for the first time the albedo and size of the very large SDO 2007 ORiq: D = 1280 ± 210 km, pr = 18.5+^-^ %. Along 
with its surface composition, this indicates that the object is Quaoar-like. 

- Our size/albedo estimations for the Eris/Dysnomia system (D = 2454 ±117 km, py = 84.5 ± 8.8 %) are compatible within error 
bars with occultation results by Sicardy et al. (2011), but further analysis is needed. 

- We find a significant correlation between albedo and diameter (more reflective objects being bigger), probably due to the fact 
that large objects can retain bright ices more easily than small objects. Stansberry et al. (2008) obtained a similar result for a 
larger sample of TNOs, but they noted that the correlation is sensitive to changes in the classification used (i.e. MPC DES; 
lElfiot et al. 20051) . 

- We also find positive correlations between albedo and perihelion distance, and between diameter and perihelion distance 
(brighter and bigger objects having larger perihelia). The first correlation has been seen before based on Spitzer measure- 
ments (Stansberry et al. 2008) and interpreted in terms of increased ice sublimation and/or increased space weathering at low 
heliocentric distances. It seems more difficult to find a physical explanation for the second correlation. 
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Table 6. Summary of correlations with diameter and geometric albedo 
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Spearman-rank correlation tests given the main correlations between the various magnitudes. *Indicates that Eris is excluded, n: number of data 
points, p: Spearman-rank correlation without data error bars. {p)V^i Spearman-rank correlation accounting for data error bars, p: significance 
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with respect to the ecliptic plane, qi perihelion, pv' V-band geometric albedo. The last line shows that rh and q are strongly correlated in our 
sample. 
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Fig. A.l. Left four panels: Typhon red band 1st visit single maps (tl to t4). Right four panels: Typhon 2nd visit single maps (t5 
to t8). tl, t3, t5, and t8 were observed with a scan angle of 110°, and t2, t4, t6, and t7 with a scan angle of 70° with respect to the 
detector array. The green circle marks the object position. 



Appendix A: Data reduction 

In this appendix, we describe the reduction process applied to the PACS data in order to obtain what we call "single" and "combined" 
maps. 

A.1. General PACS photometer data reduction 

As explained in the main text (see Sect. |3.1| ) the data reduction of the maps was performed within HIPE (version 6.0.2044). We 
modified and adapted the standard HIPE scripts for the needs of our programme. Thus, the PACS photometer data is reduced from 
Level-0 to Level-2 using a modified version of the standard mini-scanmap pipeline. A major diff'erence of our reduction chain 
compared with the standard pipeline is the application of frame selection on scan speed rather than on building block identifier (also 
known as BBID, which is a number that identifies the building block -a consistent part- of an observation). It has been widely tested 
that the application of the scan speed selection increases the number of useable frames, and finally increases the SNR of the final 
maps by 10-30% (depending on the PACS observation requests (AOR) setup and on the band used), which is especially important 
for our faint targets. 

Since our targets move slowly (a few ''/hour), and the total observation time per scanning direction in one visit is between 10-25 
min we do not correct for the proper motion. Performing this kind of correction would smear the background, and would make it 



impossible to make a correct background subtraction in the later stages (see subsection |A.2 ). 

We use a two-stage high-pass filtering procedure, as in the standard mini scan-map pipeline. In the first stage a "naive" map is 
created. This map is used in the second stage to mask the neighborhood of the source and the high flux pixels. In the second stage 
the vicinity of the source has been masked (2 times the full width at half maximum -FWHM- of the beam of the actual band), as 
well as all pixels with fluxes above 2.5x the standard deviation of the map flux values. We used the standard high-pass filter width 
parameters of 15, 20 and 35. The final maps are created by the PhotProject() task, using the default dropsize parameter. 

These first steps lead to the production of one "single" map per visit, filter, and scan direction (i.e. in total 8 maps per object 
in the red, and 4 maps in the blue/green created using the MosaicTask() in HIPE). The sampling of the single maps generated with 
these scripts are l.V /pixel, lA''/ pixel, and 2.V'/ pixel for the blue (70 jim), green (100 jim) and red (160 jim) bands respectively. 



A.2. Building final maps for pfiotometry 



We then use these single maps to generate the final "combined" maps on which the photometry will be performed. Essentially, 
to generate the combined maps, we determine the background map, subtract it from each single map, and finally co-add all the 
background- subtracted maps. A detailed description of this process is as follows: 

i) We have initially a set of 8 (red) or 4 (blue and green) single maps taken on diff'erent dates. Considering the red maps, let us 
call tl to t4 (t5 to t8) the dates corresponding to the first (second) visit. tl-t4 and t5-t8 are separated by typically 1-2 days so that the 
target motion (typically 30-50'' as described in Sect.[2]of the main text) produces a significant change in the observed sky field (see 
Figure [ATT] ). 

ii) A background map is generated using the single maps. To do this we "mask" the target in each of the 8 (4) single maps and 
co-add the map s in t he sky coordinate system. This step produces a background map with high SNR without the target (see the right 
panel in Figure [AJ] for the Typhon case at 160 micron). 



iii) The background map is subtracted from the single ma ps (Fi gure |A.l| ), producing 8 (4) single maps with background removed. 
We call these images "background-subtracted maps" (Figure [Al2|). 



iv) Finally the background- subtracted maps are co-added in the target frame (center panel of Figure [A3] ), producing the final 
combined map on which photometry is performed. 
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Fig. A.2. Left: Same as Figure 1 for the background-subtracted maps. 






Fig. A.3. Left: Typhon simple coadded map in red band. Note that in this case the background is the average (centered in the target) 
of backgrounds in maps tl-t8 of Figure [ATT] so that the bright background sources appear twice. Green circle marks the object 
position. Center: Typhon background-subtracted coadded map in red band. Right: Background map (target masked). 



Note that in step ii) before combining the images of the two epochs, we need to calculate the optimal matching offset between the 
coordinate systems of the two maps since, when we use the original World Coordinate System (WCS) information stored in the file 
headers, artifacts of bright sources remain in the background subtracted image. The best off'set position was found in the following 
way: we added a small off'set -in steps of 0.5'' in both the R.A. and DEC directions- to the appropriate original WCS parameters 
in the file headers, and performed the background subtraction. For each of these maps the standard deviation was calculated for 
pixels with a coverage value of at least 35 % of the maximum coverage of the map. The off'set values related to the lowest standard 
deviation correspond to the best matching of the two backgrounds. After this procedure the artifacts of bright sources (e.g. pairs 
of positive-negative spots) completely disappeared (or, at least, were below detection), and the noise in the background decreased 
notably. This method was performed on the 160 jim band measurements, where the eff'ect of the background is the largest. However, 
it has been checked in a series of measurements that the same off'sets can be applied to the 70 and 100 yum images as well. 

An alternate, simpler method, to obtain final maps is to co-add directly the original single maps (i.e. not background-corrected) 
in the target frame. We call this the simple coaddition method (left panel of Figure A.3 ). This method is obviously less optimal than 
the previous one in terms of SNR, but provides a useful test to demonstrate that the background subtraction is not introducing any 
spurious eff'ects. 

A more detailed description about the whole PACS data reduction process will be published in Kiss et al. {in prep.-di). 



Appendix B: Statistical aspects 

B.1. "Reseated error bar" approach 



As indicated in the main text (see Sect. 4.3 ), error bars on the fitted parameters (diameter, albedo and beaming factor) are deter- 
mined through a Monte-Carlo approach which essentially consists of randomly building synthetic datasets using uncertainties in the 
measured fluxes. However, when the measurement errors bars are too small, or similarly when the model is not completely adequate 
(both cases being indicated by a "poor model fit"), the described Monte-Carlo approach will underestimate the uncertainties on the 
solution parameters. In this case, we adopted a "rescaled error bar" approach described hereafter. 
Quantitatively, the fit quality is determined by the value of the reduced 

2 1^ (0- Mf 
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Fig. A.4. Aperture-corrected curve-of-growth obtained from Typhon background- subtracted coadded map in red band (Figure 
center). 



A.3 



where O, M and cr are the observed, modelled and error bar flux values, and v is the degree of freedom (v = N-1 for fixed-77, 
and y = N-2 for floating-//, where N is the number of thermal wavelengths available). While ~ 1 indicates a good fit^x^^^ » 
1 indicates a poor fit. The idea is therefore to empirically "rescale" (uniformly, in the lack of a better choice) all error bars cr by 

-yjx^Qd t^efore running the Monte-Carlo approach. This method leads to much more realistic error bars on the solution parameters. 

However, the range ofx^^^ for which this approach is warranted depends on the number of data points, as for few observations 
does not have to be closely equal to 1 to indicate a good fit. Specifically, the variance of the variance estimator of a distribution of 
N points picked fro m a G aussian distribution with dispersion = 1 is itself a Gaussian with a mean equal to unity and a standard 
deviation equal to ^/2/N. If we had for example N = 1000 observations, 68.2 % of the good fits have a reduced between 0.955 
and 1.045, i.e.;^^^^ is strongly constrained to 1; for A/^ = 5 (resp. 3) observations in contrast, this range is much broader, 0.368-1.632 
(resp. 0.184-1.816). This means for example that for 3 data points (PACS only), a reduced of e.g. 1.6 is acceptable. We considered 
these A/^-dependent limits as the thresholds beyond which the fits are statistically bad and the error bars need to be rescaled. Note 
finally that the rescaling approach is merely an "operative" method to avoid unrealistically low error bars on the fitted diameter, 
albedo and beaming factor, and that error bars shown in Fig. [T] and [2] of the main text are the original, not rescaled measurement 
uncertainties. 



B.2. Spearman-rank correlation with error bars 

As described in the main text the Spearman-p correlation test is distribution-free and less sensitive to outliers than most other 
methods but, like others, treats data-points as 'exact' and does not take into account their possible error bars. Any variations in 
the measured data points within their error bars, may change the correlation coefficient, however. Furthermore, each correlation 
coefficient has its own confidence interval, which depends on the number of data points and on the magnitude of the correlation. To 
account for these eff'ects we used the following three procedures: 

1 . To estimate the most probable correlation coefficient given the error bars, we generate 1000 samples of data points, building each 
synthetic dataset from its associated probability function. When errors bars are symmetric this probability function is considered 
to be Gaussian and the error bar correspond to one standard deviation. When errors bars are asymmetric, we fit a lognormal 
probability function such that the observed value corresponds to the distribution's mode, and the interval between the upper 
and lower limits of the observation corresponds to the shortest interval containing 68.2% of all possible values. The resulting 
distribution of correlation values is not Gaussian, but its Fisher-transform z = arg tanh(p) is. So, we can determine the most 
probable correlation value, (p) from our Monte-Carlo simulations, and its upper and lower limits (-\-crMc/-o'Mc), which are 
not necessarily symmetric after the reconversion using p = tanh(z) (see IPeixinho et al. 20041 for more details). It is noticeable 
(and expected) that observational error bars 'degrade' the correlation. The approximate significance p- value, or confidence level 
(CL), of p can be computed from t = p ^/(n - 2)/(l -p^), which follows a t- Student distribution with n-2 degrees of freedom. 
However, when ^ < 15, using t provides a slight overestimation of the significance. We have compared the previous calculation 
of the significance with exact values tabulated by Ramsey (1989) and computed adjustments required to obtain a more accurate 
approximation of the p's and p- values in the case of low n. 

2. The confidence level of a given p results from testing the hypothesis 'the sample is correlated' against the hypothesis 'the sample 
correlation is zero'. Knowing the confidence interval (CI) within which the correlation p value of the parent population lies may 
be more informative. For example, suppose we have p = 0.7 with a 3cr confidence level; we would conclude that the sample 
is correlated. But if the 68% confidence interval of the correlation were say [0.3,0.9], we would be unsure if the correlation 
was very strong or rather weak. To estimate the shortest (+(Tfi/-(Tfi) interval containing 68.2% of the population's p (i.e. the 
equivalent to the Gaussian Ic r interv al), we used 1000 bootstrap extractions from the data-sample computing this range as we 
did for the previous item (e.g. lEfron & Tibshirani 1 99 3 1) . 
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3. Finally, to be fully correct, one should perform the 1000 bootstraps on each one of our 1000 Monte-Carlo simulations to obtain 
the true 68% confidence interval of (p), a computationnally heavy task. Fortunately, the combination of the two eff'ects can be 
obtained by quadratically adding the stand ard deviation s of the Gaussian distributions of the Fisher-transform of the Monte- 
Carlo simulations and the bootstraps, i.e. V^mc2 + 0-52, which, after reconversion to p, will give our final best estimate for the 
68% confidence interval noted as (p)^J. 



